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Abstract 

Motivated by recent claims stating that the acceleration of the present Universe is due 
to fluctuations with wavelength larger than the Hubble radius, we present a general 
analysis of various perturbative solutions of fully inhomogeneous Einstein equations 
supplemented by a perfect fluid. The equivalence of formally different gradient expan- 
sions is demonstrated. If the barotropic index vanishes, the deceleration parameter is 
always positive semi-definite. 



e-mail address: massimo.giovannini@cern.ch 



1 Introduction 



Type la supernovae seem to suggest that the present Universe is experiencing a phase of 
accelerated expansion [1, 2] (see also [3, 4] for earlier results reporting the first evidences 
of this phenomenon). The experimental results are also consistent with the determination 
of the cosmological parameters coming from CMB physics [5] . These results seem to imply 
that up to 70 % of the present energetic content of the Universe is "dark" . If the dark 
energy component is parametrized with a perfect fluid characterized by a barotropic index 
w, the experimental determinations reported in Refs. [1, 2] suggest that for a flat Universe 
with Q tot = 1, ft mattcr = 0.29±8;{jjj and w = -1.02±{J;$. 

The simplest model of dark energy is, in some sense, a time-independent cosmologi- 
cal constant whose associated w is exactly —1. Fluid models of dark energy may lead to 
barotropic indices in the range —1 < w < —1/3. It is equally plausible that the dark-energy 
component is described by a scalar degree of freedom (i.e. the quintessence field) whose 
potential becomes dominant around the present time [6] (see also [7] ) . Such a degree of free- 
dom may also be a pseudo-scalar [8]. Quintessence models of dark energy imply generically 
that w > — 1. A rather economical class of models is the one where the quintessence field 
(dominant today) and the inflaton field (dominant in the far past history of the Universe) 
are identified in a single (scalar) degree of freedom. These are the so-called quintessential 
inflationary models [9] whose salient features may lead necessarily to a stochastic back- 
ground of gravitational radiation for typical frequencies larger than 0.1 KHz [10, 11, 12]. 
Useful surveys on the status of fluid and quintessential models of dark energy can be found, 
respectively, in Refs. [13] and [14]. 

Experimental data cannot rule out values w < — 1: this occurrence stimulated the 
investigation of fluid models of dark energy where the cosmological energy density is future 
increasing, rather than future decreasing, as in the case w > — 1. Provided w < — 1 all 
the time in the future, a "big rip" can be expected [15], i.e. an infinite expansion of the 
Universe in a finite amount of time [16, 17]. Along this perspective, recently Barrow has 
discussed the interesting case of "sudden" singularities that may arise in expanding FRW 
Universes even if the dominant energy condition is not violated [19, 20]. The possibility 
of anisotropic sudden singularities has been also discussed recently by Barrow and Tsagas 
[21]. We will briefly comment about the possible relevance of the considerations developed 
here in the context of the physics of sudden singularity. 

The aim of the present paper is to understand wether an inhomogeneous Universe may 
sustain a phase of accelerated expansion even if the barotropic index of the sources is pos- 
itive semi-definite, i.e. w > 0. There are various reasons to discuss this problem. It is, 
in fact, plausible that inhomogeneities (with typical wavelengths larger than the Hubble 
radius) may be generated in the early Universe provided the duration of the inflationary 
phase exceeds, say, 65 e- folds. After equality and before decoupling, (when inflationary 
perturbations are imprinted, via the Sachs- Wolfe effect, on the temperature anisotropies) 
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the metric fluctuations are perturbative. Still, one could wonder if higher-order effects in 
the amplitude of the fluctuations may change this perspective and make the present ampli- 
tude of inflationary generated fluctuations so large to affect, for instance, the deceleration 
parameter. 

There may be a way out in this technical impasse. The idea is to forget about homo- 
geneous backgrounds (supplemented by tiny inhomogeneous fluctuations) and to consider, 
instead, some type of gradient expansion. This theoretical tool has a long history that 
can be traced back to the seminal contributions of Lisfhitz, Belinskii, Khalatnikov and 
collaborators (see [22, 23] and references therein). This approach allows to compute the 
phenomenologically relevant quantities to a given order in the spatial gradients. There have 
been recently claims [24] suggesting that indeed a matter-dominated Universe (without any 
vorticity) may be accelerating solely thanks to the presence of super-Hubble inhomogeneities 
2 . In [25, 26, 27, 28, 29, 30, 31] various independent criticisms of this proposal have been 
presented. In [25] it was argued that fluctuations with typical wavelength larger than the 
Hubble radius can certainly affect the spatial curvature but can never accelerate the Uni- 
verse. Similar conclusion, through a different chain of arguments, has been derived in [26]. 
In [27] demonstrated that super-horizon fluctuations do not produce an accelerated expan- 
sion. Their approach has been to take an exact solution of Einstein equations 3 (describing 
an underdense region). This strategy allowed to obtain a specific (computable) model of 
super-horizon fluctuations falling in the class discussed in [24]. The example showed that 
no acceleration takes place in agreement with general theorems implying that the Universe 
cannot accelerate unless the strong energy condition is violated. In [27] an explicit analysis 
of the variance of the deceleration parameter has been also performed in the case of single- 
field inflationary models. In [29] the perspective has been to analyse a form of averaged 
Einstein equations by assuming the correctness of some of the conclusions of [24]. The 
conclusion of this study has been that the model of [24] is likely to be ruled out by present 
observations. 

The analysis of Ref. [30] dealt also with the effect of small-scale fluctuations (i.e. with 
wavelengths smaller than the Hubble radius). The conclusion of the analysis has been that 
sub-horizon perturbations are not a viable candidate for the explanation of the present 
acceleration of the Universe. 

In [31], an approach for the analysis of fully inhomogeneous Einstein equations has been 
proposed and applied to the case of Ref. [24]. The calculation of [24] has been then repeated 
and it has been discovered that, within the accuracy of the solution, the late-time behaviour 
of the spatial gradients in a matter dominated Universe can never make the deceleration 
parameter negative. Furthermore, by enforcing the validity of the perturbative expansion it 

2 To stress the fact that the Universe is accelerated thanks to super-horizon fluctuations, the authors of 
[24] called their proposal Super Hubble Cold Dark Matter (SHCDM) to be contrasted with the usual ACDM. 

3 Recently, a similar perspective has been adopted in the analysis of a Lemaitre-Tolman-Bondi Universe 
by Alnes, Amarzguioui and Gron [28] (see also [18] ). 
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has been shown that, at later times the deceleration parameter vanishes, i.e. the expected 
result is the Universe is dominated by gradients. 

The purpose of the present paper is to extend and generalize the analysis presented 
in [31] along different lines. The first generalization consist in discussing the gradient 
expansion without the restriction of a dust-dominated Universe but for a generic barotropic 
index. Furthermore, in [31] the analysis has been performed by only keeping only one scalar 
degree of freedom of the fully inhomogeneous metric. Here, as announced [31] all the degrees 
of freedom of the inhomogeneous metric will be analyzed simultaneously. The related aim 
of the present analysis is to stress the mutual connections of apparently diverse gradient 
expansions. 

The plan of the present paper is the following. In Section 2 the fully inhomogeneous 
Einstein equations will be introduced in terms of the extrinsic and intrinsic curvature ten- 
sors. Section 3 deals with the interesting case of anti-newtonian solutions that could be an 
appropriate seed solution for a gradient expansion aiming at the description of super-horizon 
fluctuations. In Section 4 the gradient expansion will be derived. Solutions accurate up to 
two gradients will be presented. In Section 5 the gradient expansion will be applied to a 
situation where the parametrization of the seed metric is apparently different. Section 6 
contains some more phenomenological considerations dealing with the late-time behaviour 
of a Universe dominated by super-Hubble fluctuations. Finally, Section 7 contains some 
concluding remarks. 



2 Inhomogeneous Einstein equations 

Consider a four-dimensional line element of the form 

ds 2 = dt 2 - jij{t, i)diW; (2.1) 

the symmetric rank-two tensor 7^ in three dimensions contains 6 independent degrees of 
freedom. The components of the extrinsic curvature are 

In the following, the compact notation TrK 2 = K\ Kj will be used. The three-dimensional 
Ricci tensor becomes instead 

_ q r m — f> v m j_r m r i _ r m (0 *x\ 

1 13— u m i - ij U J L mi 1 L ij L mi L jm L iti \^- ) 

where 

r™ = ^ me (-daij + d ]lu + d ilje ). (2.4) 

Equations (2.2) and (2.3) allow to write the Einstein equations in a fully inhomogeneous 
form. More specifically, assuming that the energy-momentum tensor is a perfect relativistic 
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fluid 

z£ = (p + pK«"-p^. ( 2 - 5 ) 

the Hamiltonian and momentum constraints are, respectively, 

K 2 - TrK 2 + r = 16irG[(p + p)u u° - p] , (2.6) 
ViK- V k K? = 8irG Ui u°(p + p). (2.7) 

The (ij) components of Einstein equations lead instead to 

{ kf - K K{ - KSf) + \81{K 2 + TrK 2 ) - (r\ - \r%) 

= -8irG[(p + p)u iU j +p8j], (2.8) 

where the overdot denotes a partial derivation with respect to t while Vj denotes the 
covariant derivative defined in terms of jij and of Eq. (2.4). A trivial remark is that, in 
Eqs. (2.6), (2.7) and (2.8), the indices are raised and lowered using directly jij(t,x). 

By combining the previous set of equations the following relation can be easily deduced 



qTrK 2 = 8vrG 

where 



(p + p)u u° + 



p- p 



(2.9) 



= -l+n§2. ( 2 - 10 ) 

is the inhomogeneous generalization of the deceleration parameter. In fact, in the homo- 
geneous and isotropic limit, 7^ = a 2 (t)5ij, K\ = —Hdj and, as expected, q(t) — > —da/a 2 . 
Recalling the definition of TrK 2 (given after Eq. (2.2)) it is rather easy to show that 

K 2 

TrK 2 > > 0, (2.11) 

3 

where the sign of equality (in the first relation) is reached, again, in the isotropic limit. 
Since 7^ is always positive semi-definite, it is also clear that 

u u° = 1 + j ij UiUj > 1, (2.12) 

that follows from the condition g^u^u u — 1. From Eq. (2.9) it follows that q(t,x) is 
always positive semi-definite if (p + 3p) > 0. As correctly pointed out in [27] and [29] 
that this statement assumes the absence of rotational fluctuations. This is rather plausible 
since super-horizon vector modes are very unlikely to be generated either during inflation or 
during a phase where the Universe expands. Concerning the possible role of vector modes 
in conventional and unconventional inflationary models see, for instance, [32, 33, 34, 35]. 
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Equations (2.6), (2.8) and (2.7) must be supplemented by the explicit form of the co- 
variant conservation equations: 



-^[V7CP + P> Q u l ] - -±=d k { V7[(P + P)nV + pj ki }} - 2K\u°i/{p + p) 



Id 1 

-y^iv^KP + p)uou° -p}}~ ~^di{^(p + p)u u 1 } 



(2.13) 
(2.14) 



-K e k [(p + p)u k u e +pS$] =0, 
where 7 = det(7y). 

It is useful to recall, from the Bianchi identities, that the intrinsic curvature tensor and 
its trace satisfy the following identity 



V,-r: 



1. 



3 i 



-V,T. 



(2.15) 



Note, finally, that combining Eq. (2.6) with the trace of Eq. (2.8) the following equation 
is obtained: 

TrK 2 + K 2 + r-2K = 8wG(p-Sp). (2.16) 
Equation (2.16) allows to re- write Eqs. (2.6), (2.8) and (2.7) as 



K - TrK = 8itG 



( _l ^ , P-P 
(p + p)u u H — 



1 d 



8vrG 



ViK - V k K k = 8ttG( P + p)ui u°. 



(p + p)UiU> + ^-y^^ 



(2.17) 

(2.18) 
(2.19) 



3 Anti-newtonian solutions 

Consider the situation where the scalar curvarure r and the velocities are all much smaller 
than K 2 and K. This situation describes the occurrence of a "anti-newtonian" regime and 
has been previously investigated by Tomita [36, 37] as well as by Deruelle and collaborators 
[38, 39, 40] interested in various aspects of this approximation. 

If K 2 and K are both leading in comparison with the curvature and velocity contribu- 
tions, Eqs. (2.17), (2.18) and (2.13) can be written, as 

K 2 

K TvQ 2 = At:G(3p + p), (3.1) 

~(^K) = 12*G(p - p), (3.2) 

^|0/7« = 0, (3.3) 
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where we defined 

Kl = ^+Ql (3.4) 

with Q = Q\ = 0. From Eq. (2.14), covariant conservation of the energy-momentum tensor 
implies 

P + ^(P + P) = 0- (3-5) 

From Eq. (3.3) and (3.5) it follows 

/ 9 = / 9 (x)7~^, (3-6) 
Ql = ^ (3.7) 

where A^ = and A^ only depends on space (and not on time); note that we used the fact 
that p = wp. 

Inserting then Eq. (3.2) into Eq. (3.1) to eliminate the energy density we obtain: 

*-H±i*» + f(„ -1)^=0, (3.8) 
where, as usual, TrA 2 = A^A*-. A new variable can now be defined, namely, 

w+l 

M = -f—; (3.9) 
recalling that 7 = —27 If, Eq. (3.8) becomes, in terms of M: 

On r» w — 3 

M = -(w - l)TrA M~. (3.10) 

8 

Integrating once with respect to the cosmic time, the following relation can be obtained: 

M 2 = hw + l) 2 TT\ 2 M 2J ^ +N(x). (3.11) 
8 

Equation (3.11) can be solved (either analytically or numerically) for a given value of the 
barotropic index. In particular, in the case w = 0, recalling Eq. (3.9) we can write: 

3 TrA 2 

^=N(x)[t-t 1 (x)} 2 -^ J ^- y (3.12) 

where a further integration constant appear. Either N(x) or t\{x) can be fixed by exploit- 
ing the remaining gauge freedom of the synchronous system. Note that the momentum 
constraint of Eq. (2.19) implies 



U°Ui — 



12irGp(w + 1 

that is of higher order in the gradients, as expected 



3 k t k 

- ^-^=(r^Aj - r im A™) 



(3.13) 
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The deceleration parameter defined in Eq. (2.10) can be written, using Eqs. (2.10) and 
(3.7) as 

_ , 1 f2jK 2 - 3TrA 2 \ , t . 



i.e., using Eq. (3.12) 



_ 1 / 8JV(x) 2 [t-t 1 (x)] 2 -3'ftA 2 \ 
? ^ " 4UA r (^ 2 [t-ti(^)] 2 + 3TrA 2 J' 1 j 



whose large time limit is 1/2. Note that by looking superficially at Eq. (3.15) one would 
be tempted to conclude that q might be negative. This is not correct since the numerator 
in of the expression in Eq. (3.15) is exactly proportional to ^ > 0. 

Now one could choose, as seed metric, the solution found with this method. Then 
compute the spatial curvature and obtain the following order in the gradient expansion and 
so on [38, 40]. However, the solution of Eq. (3.11) cannot be inverted analytically for a 
generic barotropic index w. 



4 Quasi-isotropic solution and gradient expansion 

The solutions of Eqs. (2.6)-(2.7) can be classified, according to their degree of isotropy, in 
quasi-isotropic solutions and fully anisotropic solutions. For instance, close to the initial 
(big-bang) singularity the solution of Eqs. (2.6)-(2.7) is, in general, fully anisotropic. In the 
opposite limit, i.e. far from the initial singularity, the possibility of quasi-isotropic solutions 
becomes more relevant. Quasi-isotropic solutions exist, indeed, only in the presence of 
matter [41, 42] (see also [43, 44, 45]). 

Let us now look for solutions of the previous system of equations in the form of a gradient 
expansion. In other words we shall be considering 7^ written in the form 

llk = a 2 {t)[a lk {x) + (3 ik (t,x)], ^ = ~ kj (t,x)], (4.1) 

where (3(x, t) is considered to be the first-order correction in the spatial gradient expansion. 
Note that from Eq. (4.1) 7^7^ = 6j + 0{(3 2 ). The logic is now very simple since Einstein 
equations will determine the specific form of once the specific form of «y is known. 

Using Eq. (4.1) into Eqs. (2.2) we obtain 

Kl = -(H5{ + ^y K = -(3H+±0), TrK 2 = 3H 2 + Hp, (4.2) 

where, with obvious notations, H = a/a. 
From Eq. (2.7) it also follows 

V fc # - Vi/3 = IGttGu; u° (p + p) . (4.3) 



7 



The explicit form of the momentum constraint suggests to look for the solution in a separable 
form, namely, (3f(t,x) = g(t)Bj(x). Thus Eq. (4.3) becomes 

g{V k Bf - V l B) = 167rGu iU °(p + p). (4.4) 

Using this parametrization and solving the constraint for m, Einstein equations to second 
order in the gradient expansion reduce then to the following equation: 

2V-' w — 1 

(g + 3Hg)B\ + HgBSj + ^ = ^jCq + 2Hg)B5{. (4.5) 
In Eq. (4.5) the spatial curvature tensor has been parametrized as 

V j 
a 

Recalling that 



r\ = -L. (4.6) 



H = H a- SJ ^ 1 , H = -^±^H\ (4.7) 



the solution for Eq. (4.5) can be written as 



& = t (pi - 5 + 6^-3^ 2 ■ x 

H$(3w + l)(3w + 5) v 1 4(9w; + 5) ' ' ' 

with g(t) simply given by 

g (t)=a 3w+1 . (4.9) 

Note that, in Eq. (4.8), Hq = 2/[3(w + 1) to]- Eqs. (4.8) agree with the expression obtained 
in the case w = 1/3 in Ref. [22] (where the notations are such that t = 1). Equation (4.8) 
can be also inverted, i.e. Vf can be easily expressed in terms of B\ and B: 

tt2 

V\ = — ^[B5j(6w + 5 - 3w 2 ) +Bl(3w + 5)(3w + 1)] (4.10) 
Using Eq. (2.15) the peculiar velocity field and the energy density can also be written 

as 

3 / w 



vPllj = — - 



^ja 3w+1 HdiB(x), 



3H% 



8vrGpV3w + 5 
'a-^+V - ^±±B(x)a 



r 8ttG 

Let us therefore rewrite the solution in terms of 7^, i.e. 

aij (x)+a 3w+1 B ijy 



lij = a (t)[ctij(x) + @ij(x,t)] = a (t) 
Concerning this solution a few comments are in order: 



(4.11) 



(4.12) 
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• if w > —1/3, Pij becomes large as a — > oo (note that if w = —1/3, a 3w+1 is constant); 

• if w < —1/3, flij vanishes as a — ► oo; 

• if w < — 1, not only the gradients become sub-leading but the energy density also 
increases as a — > oo. 

• to the following order in the perturbative expansion the time-dependence is easy to 
show: 7^ ~ a 2 (t)[aij + a 3w+1 Bij + a 2( - 3w+1 ^£ij] and so on for even higher order terms; 
clearly the calculation of the curvature tensors will now be just a bit more cumbersome. 

The first two comments are rather elementary. The first comment simply expresses the fact 
that if w > — 1/3 the gradients decay close to the initial big-bang singularity (but not far 
from it!). The second comment simply justifies why during a phase dominated by an effective 
cosmological constant (or by a fluid violating the strong energy condition) the gradients are 
washed out. The third comment may have some interesting implications for the study of 
big rip singularities or for the more general case of sudden [19, 20, 21] singularities. It would 
be for instance interesting to investigate if, in general, sudden singularities will exhibit or 
not some type of BKL oscillations that are known to be present in the case of the initial 
big-bang singularity. 



5 Equivalent forms of gradient expansion 



Consider now the following parametrization of the perturbed metric: 



-2¥(t,2)_2 



a (t)[5ij + mj(x,t)], 



(5.1) 



where p\ = 0. This parametrization has been used in [24] and in [31]. Equation (5.1) may 
seem, superficially, different from the one of Eq. (4.1) since, on one hand, the analog of 
a.{j has a very specific tensor structure coinciding with the three-dimensional Kroeneker 
symbol; on the other hand VP is allowed to depend both on space and time. 

Equation (5.1) 



In this case, the Hamiltonian constraint of Eq. (2.6) leads to 



„2* 



6(1r-tf) 2 + — - 



4V 2 ^ - 2(W) 2 ] = WttG 



P+{P + P)— o"(« 



u 2 ) 



(5.2) 



(5.3) 



where u 2 = UiUj5 1 ^ and it 2 = UiUj/j, 1 ^. From the (ij) components of Einstein equations we 
get, after linear combination, 



2* 



3(tf -H)- 3(tf - H) = 4ttG( P + 3p) + 8ttG( P + p)^-{u 



2 u 2 ) 
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= 2* 



i li+3Hf i i = -2— 
„2* r 



1 



1 



+16irG(p + p)- 



(UiU 3 — UiU 3 ) 



u 



u 2 



(5.4) 



where V 2 ^ = S^didjty and where (V^) 2 = 6 lJ d^dj^ . Finally, the momentum constraint, 
to first-order in the gradient expansion, is 



+ -dktf = 4?rG(p + p)u iU °. 



(5.5) 



A more tractable form of the system can be obtained by eliminating the energy density; 
hence the following pair of equations emerges: 



■i S(w + 1, T 



3a 2 



V 2 ^ 



0, 



A? + 3 #Ai = — 2 e 



2* 



- ^(V 2 ^)^ + diVd** - ,(W) 2 ^ 



The solution of Eqs. (5.6) and (5.7) can be written as: 



77 2 



v 2 /-^(v/) 2 



where 



" H 2 



ji(w) 



a 3»+l e 2/ 



w - \(^ 2 m + dif&f - \{vf) 2 5{ 



(5.6) 
(5.7) 

(5.8) 
(5.9) 

(5.10) 



3(9w + 5)' JZV ~ y (3w + l)(3u; + 5) 
and where f(x) is an arbitrary function depending upon space but not upon time. In this 
function is encoded the information on the specific kind of super-horizon fluctuations. The 
velocity fields have been consistently neglected since they are of higher order in the gradient 
expansion. This aspect can be appreciated by inserting the solution (5.8) into Eq. (5.5): 
di^ only receives contribution from the second term in Eq. (5.8) (i.e. the one proportional 
to ji{w)). Thus ui will contain three gradients and will be negligible at this order. Notice, 
however, that the contribution of the velocity field must be necessarily taken into account 
when going to orders higher than the second in the gradient expansion. 

Since Eq. (5.8) and (5.9) are derived under the approximation that terms with more 
than two gradients are neglected, the terms \& 2 has to be negligible with respect to ^, i.e. 



2 1* | > S(w + 1)|^| 



(5.11) 



This aspect can be appreciated by computing ^ 2 from Eq. (5.8): the only contribution 
to VP comes from the term proportional to ji(w) in Eq. (5.8) and this has four spatial 
gradients. A relevant consequence of this basic observation is that when computing TtK 2 
from Eq. (5.2) the correct result, within the approximations made so far, is 



TrK = H 2 - 2H$> > 



(5.12) 
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and not TrK 2 = H 2 — 2H^ + VP 2 . Notice that the inequality appearing in Eq. (5.12) is a 
direct consequence of the inequalities reported in Eq. (2.11): TtK 2 is always (i.e. at any 
order) positive semi-definite. This simple observation has a simple consequence, i.e. that 



* < 



H 



For future convenience, let us rewrite the solution of Eq. (5.8) as 

tf(x,t) = f(x) +a 3w+1 ^ (x) 
where can be read-off from Eq. (5.8) 



(5.13) 



(5.14) 



(5.15) 



From Eq. (5.14), by taking the first time derivative, it follows easily that \l/ = (3w + 
l)a 3w+1 H^>o, and this implies: inserting the obtained expression for 4f into Eq. (5.13), the 
inequality puts a condition on ^q, i.e. 



,3w+l 



2(3w + 1)' 



that translates, in the case w = into 



(5.16) 



(5.17) 



We are now in condition to clarify the relation between the perturbative expansion 
discussed in the present Section and the general quasi-isotropic gradient expansion derived 
in Eqs. (4.8) and (4.12). The conclusion will be that, if the positivity of TrK 2 is enforced, 
the two expansions are exactly equivalent, in the sense that the expansion discussed in the 
present Section is just as sub-case of the general quasi-isotropic expansion. 

For this purpose let us then consider Eq. (4.12) and let use the freedom to specify ctij 
in Eq. (4.12) and let us choose 

a ij (x) = e- 2 ^S ij . (5.18) 

Then, the spatial curvature tensor can be immediately computed using Eq. (5.18) and Eq. 
(2.3). Thus, recalling Eq. (4.6), V\ can be easily computed: 



=2/ 



V z f8l+d i &>f-{Vf) 2 5i+d i fd>f 



(5.19) 



Consequently, from Eq. (4.8), the tensor Bij will be 
4e 2 / 



Bl = - 



H$(3w + l)(3w + 5) 



3w(w + l) v2/ 3w 2 + 12w + 5 (v/)2 



(9w + 5) 



2(9w + 5) 



%+di&>f+dif&>fy 

(5.20) 
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This expression is still difficult to compare with the parametrization employed in the present 
section since B\ is not traceless. Let us then separate the traceless contribution by writing 

jij(x,t) as 

7lj (x,t) = a 2 (t)e- 2 *[^- +fMj{x,t)] (5.21) 

where p,\ = 0. The traceless tensor fi^ is essentially given by the traceless part of = 
a (3w+i)g_ w \ iere fy. i s given by Eq. (5.20) . Therefore, we can also write the solution in 
terms of VP: 



(5.22) 



In the specific case w = the previous expression becomes 

* = f(x) - I In { [l - ^ e 2 /a[v 2 / - i(V/) 2 ] }. (5.23) 

Let us now compare Eqs. (5.22) and (5.23) with Eq. (5.8). For for sake of simplicity let us 
consider the case w = 0. Therefore, in this case, * of Eq. (5.23) can be written as 

* = /(*) -1 In [l-2a* ], (5.24) 

where ^o(x) is the same one defined in Eq. (5.15) and it includes the contribution of the 
spatial gradients. Recalling now Eq. (5.17) we clearly see that the argument of the logarithm 
appearing in Eq. (5.24) is only compatible with TtK 2 being positive semidefinite, if 2a^o < 
1. This implies that, as anticipated, VP = * within the accuracy of the approximation. 

The condition coming from TrK 2 > is necessary for the consistency of the approxi- 
mation of the present section. It is however not sufficient. This means that an even more 
restrictive condition on V&o stems from Eq. (5.11). Indeed, using the parametrization (5.14) 
it follows, in the case w = 0, that 

a|* | < g- (5-25) 

This inequality follows directly from Eq. (5.11) by recalling, from Eq. (5.14), that, for 
w = 0, VP = / + a^o- Thus \t = (H + H 2 )\^o\a; using the last expression in Eq. (5.11) and 
recalling that 2H = —3H 2 (in the case of dusty matter) we Eq. (5.25) is readily obtained. 
As anticipated the condition (5.25) is more restrictive than the one derived in (5.17) always 
in the w = case. In fact, Eq. (5.25) implies —1/3 < a^o < 1/3 while Eq. (5.17) implies 
only that aH < 1/2. 

The last step of our discussion is to show that not only * = * but also that fiij = fiij. 
This is immediate recalling that jlij is constructed from the traceless part of fyj whose 
tensor structure in given, in terms of the intrinsic curvature tensor, by Eq. (4.8). Bearing 
now in mind Eq. (5.18), it is immediate to show that fiij equals ^ given in Eq. (5.9). 
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6 Deceleration parameter (s) 



In this section we will be concerned mainly with the calculation of the deceleration parame- 
ter. In particular we want to show that the deceleration parameter for a matter-dominated 
Universe can never become negative. In a complementary perspective one could say that 
if the deceleration parameter becomes negative, within the gradient expansion discussed in 
the present paper, the conditions of validity of the perturbative expansion must be violated. 

Let us therefore start from the general definition of the deceleration parameter valid in 
the inhomogeneous case and reported in Eq. (2.10). Then let us consider the parametriza- 
tion of Eq. (5.14) and in the case w = 0. Inserting Eq. (5.14) into Eq. (2.10) we obtain 

1 , 3/2-q^ /2 3a^ + 1 

q(x,t) = —l + — = J — = — =— . (6.1) 

Hy ' (l-2a* ) 2(l-2a* ) 

In order to derive this result, we inserted Eq. (5.14) into K, K and TrK 2 appearing in 
the general definition (2.10). Then we used the relation 2H = —3H 2 (valid in the dust- 
dominated case). Finally, we consistently neglected the terms \E' 2 that, as explained in the 
previous Section, are of higher order since they contain 4 gradients. This means, in practice, 
that TrK 2 = H 2 (l — 2a^o). Concerning Eq. (6.1) the following comments are in order: 

• as derived in the previous Section the validity of the perturbative expansion implies 
in the case w = that aj^ol < 1/3 (see Eq. (5.25) and comments therein); 

• therefore, it is clear from the expression of reported in the second equality of Eq. 
(2.10) that q is always positive semi-definite for —1/3 < a^o < 1/3; 

• furthermore, when a^>o —> —1/3 (corresponding to the maximally underdense Uni- 
verse allowed by the validity of the perturbative expansion), then q — ► 0. 

The last occurrence corresponds to a gradient-dominated Universe. This result was already 
derived in [31]. Therefore, the validity of the perturbative expansion forbids a negative 
deceleration parameter and, therefore, the acceleration of an Universe dominated by dusty 
matter (i.e. w = 0). 

Notice also that Eq. (5.25) implies an upper limit on a(t), i.e. 

1 , x 

a < (6.2 

" 3|tf | 

If sufficiently small values of are allowed, a can be rather large. To support their 
argument, the authors of Ref. [24] took (in their figure) the smallest value of to be 
~ —1/4. This simply means, according to Eq. (6.2) that a < |. In light of the discussion 
of the present paper the soundness of this results stems directly from the equivalence of the 
gradient expansions discussed in Sections 4 and 5. In fact, it is clear that the quasi- isotropic 
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gradient expansion of Section 4 is suitable in the limit a — > 0, i.e. in the limit where it was 
correctly applied in the early sixties. 

From the results of Section 4 and 5, moreover, it is clear that there are cases when the 
gradient expansion may be suitable to describe the large times behaviour. But these cases 
are precisely the ones for which w < — 1/3. In this class of fluid models, as remarked in 
Section 4, the contribution of the gradients becomes progressively less important as the 
Universe expands. Unfortunately, the cases w < —1/3 correspond to an Universe that is 
not dust-dominated. 

Finally, in order to complete our discussion, let us see what happens if we include terms 
of higher order in the gradient expansion and take the limit for a — > oo. In this case we 
have that TrK 2 = (H — 4>) 2 . Then using the parametrization of Eq. (5.14) into Eq. (2.10) 
we do find that the deceleration parameter is given by 

3/2-o*o/2 

9(M) = " 1+ (1-a^o) 2 ' (6 - 3) 
The conclusion of Ref. [24] is that if we extend the validity of Eq. (6.3) into the far 
future, the deviation of the deceleration parameter from q = 1/2 becomes larger and larger 
approaching the asymptotic value q = — 1. On the basis of the gradient expansion discussed 
in the present paper such a chain of arguments is arbitrary. 



7 Concluding remarks 



In the present paper the gradient expansion has been discussed for a generic perfect fluid 
with barotropic index w. It has been shown that formally different realizations of the 
gradient expansion are in fact equivalent if the perturbative expansion is correctly handled. 
The aim of the present analysis was to understand recent claims suggesting that perfect 
fluids with w > can act as effective dark-energy candidates. The results of the present 
study imply that if < w < 1 the deceleration parameter is always positive semi-definite. In 
the present analysis we kept terms up to two gradients in the inhomogeneous expansion. At 
this order, the negativity of the deceleration parameter for matter sources with < w < 1 
is just a signal of the breakdown of the perturbative expansion. The techniques discussed 
in the present paper may also be relevant in order to study a possible inhomogeneous (and 
anisotropic ) approach of big rips and sudden singularities. We leave this problem for future 
investigations. 



Note added. After completing the present paper we became aware of astro-ph/0506534 
by E. Kolb, S. Matarrese and A. Riotto. These authors study solutions of the gradient 
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expansion discussed in Section 5 of the present paper. Their results are less general since 
they only consider a dust-dominated Universe. We checked that our solutions coincide with 
their solutions in the limit w = (our f(x) becomes, in their notations, 5/3ip(x)). They 
also claim that the perturbative expansion they use (i.e. the one we discussed in Section 5) 
sensibly differs from the gradient expansion we discussed in Section 4 of the present paper. 
Our results show the opposite. 
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